A novel exact cosmological solution of Einstein equations 
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We present a novel homogeneous and geometrically flat exact solution of Einstein's 
General Relativity equations for an ideal fluid. The solution, which describes an 
expanding/contracting hypercylinder, fits well with the observational pillars upon 
which rely the standard FLRW cosmology and, furthermore, it can naturally solve 
some of its most outstanding problems. 

The Friedmann-Lemaitre-Robertson-Walker(FLRW) exact solution of Einstein's General 
relativity equations for an ideal fluid describes an expanding/ contracting homogeneous and 
isotropic universe and, hence, it implements the so called Cosmological Principle. This 
cosmological model, regarded in the literature as the standard cosmology, relies on three solid 
observational pillars: a) the relation between redshift and distance of nearby astrophysical 
objects, b) the extremely isotropic microwave blackbody radiation that permeates the whole 
sky and c) the collected data on primordial densities of light nuclei. The model, whose metric 
element is given by ds^ = dt^ — o?{t) (j^^^ + r'^dD?^, considers three possibilities: a) = 
for a spatially flat euclidean universe; b) A; > for a closed (hyper)spherical universe and c) 
A; < for an open hyperbolical universe. The available observational precision data favors 
the first option, a geometrically fiat universe. 

Nonetheless, it is widely accepted that the standard FLRW cosmology is an incomplete 
model. Some unsolved theoretical issues (the flatness problem, the horizon problem,...) and 
disagreements between its predictions and the observational data (the baryon density needed 
to fix the data on primordial densities of light nuclei barely accounts for 5% of the critical 
mass density needed for a geometrically fiat universe) have forced cosmologists to incorporate 
into the standard FLRW model new speculative ingredients. Namely, a very early stage of 
cosmological inflation is widely assumed to have preceded the standard evolution of the 
universe, even though there does not exist yet neither theoretical nor direct observational 
confirmation of the mechanism that could have driven it. Furthermore, it is necessary to 
postulate that some mysterious dark energy or cosmological constant of unknown nature 
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accounts for almost 70% of the energy content of the universe, while another 25% is due 
to some kind of cold dark matter. This phenomenological extension of the standard FLRW 
flat cosmology is known in the literature as concordance ACDM cosmology, for it largely 
reproduces the precision data collected during the last twenty years [HE]. Yet, some authors 
have noticed intriguing disagreements [3] and also alternative cosmological scenarios has been 
proposed [1]. 

Here we present a novel spatially homogeneous and geometrically fiat exact solution of 
Einstein's General Relativity equations and explore the possibility that it could describe 
the patch of the universe observable to us. This novel solution, which describes an expand- 
ing/contracting hypercylinder, naturally solves the flatness problem of standard FLRW cos- 
mology. In addition, it can fit the observational data upon which relies the standard FLRW 
cosmology without any need to postulate any kind of mysterious dark energy. 

We consider a 4-dim Lorentzian manifold whose metric in a cartesian grid x = 
{x^ , , x"^ , x^) we parameterize as 



goo{x) = 1, goi{x) = gio{x) = 0, gij{x) 



(1) 



where t = x° is cosmological time, {x*} 



i=l,2,3 



are comoving coordinates over the 3-dim 



spatial sheet, r = [dijX^x^)^^"^ is a comoving radial coordinate, f2j 



6ij^ are 



normalized spatial coordinates and tq fixes a reference comoving scale. This metric describes 
and expanding/contracting hypercylinder with comoving radius tq (see Fig.l). This can be 
immediately noticed by writing the invariant metric element as 



ds^ = g^.^clx'^dx" = dt^ - [a^{t) {Q,dx') {Qjdx^) + b^{t) {6ij - QiQj) dx'dx^] 
= de - rl [a\t)df^ + h\t) {de^ + sin^ed(f)^)] , 

where f = ln{r /tq) and ^, are comoving spherical polar coordinates. In order to obtain the 
last expression we write dx* = d{r fi*) = r dfi* + dr and then notice the normalization con- 
strain dijVf^Vt^ = 1, which implies 6ijQ'^dQ^ = 0. Hence, the dimesionless scale factors a{t) 
describes the expansion/contraction of the hypercylinder along the radial (axial) direction, 
while b{t) describes the expansion/contraction of its physical radius. 

After some straightforward calculations we can find the components of the Ricci tensor 
associated to the given manifold. 
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Roo{x) = Ro{t), Roi{x) = Rio{x) = 0, Rij{x) = ^ + {Sij - 0,%)) ,(3) 

where Ro{t) = - (|| + 2||), = a{t)d{t) + 2a{t)a{t)'^^ and Reit) = b{tjb{t) + 

b{ty + ^b{t)b{t) + r^^ are functions of cosmic time only. The corresponding Ricci scalar 
R — RfiuQ^'' is constant over the spatial sheet, 

^^""^ ~ ^ [a{t) + ^b{t) + \it) bit) ^ b{tr ^ b{ty ) • 

The components of the Einstein tensor G^,^ = R^^, — ^Rg^u are, therefore, given by: 



Goo(x) = Go(t), Goi(x) = Gio(x) = 0, Gij(x) = ^ (G'A(t)Jl^^- + Geit) (S^j - Q,^)) ,(5) 
where Go{t) ^ 2|lig + gj^ + ^ - (2|| + + ^) a\t) and ^ 

\^a(i) ^ fe(t) a(t) y " VV- 

For an ideal rclativistic fluid the energy-momentum tensor is given by T^yi^x) — 
\e{x) ^ P{xy\Up{x)ui,{x) — P{x)gf^i,{x), where e{x) and P{x) are, respectively, the energy 
density and pressure of the fluid in its proper local frame of reference and is the velocity 
four- vector. If we assume that the fluid is at rest in the comoving grid, then — (1,0,0,0) 
and the energy-momentum tensor takes the form. 



Tooix) = e{x), Toiix) = T^oix) = 0, Tij{x) = P(a;) {a\t)fliflj + b\t) {Sij - QiQ 



(6) 



Therefore, Einstein equations Gfj^i, + Ag^i, = SnGT^iy, where for the sake of generality we 
have included a non zero cosmological constant A, can be exactly solved if both the energy 
density and the pressure of the fluid are homogeneous over the spatial sheet and evolve in 
time according to the set of equations 



e{t) 



Pit) 



1 



SttG 



a{t) b{t) b{t)^ r, 



+ 



+ 



a{t)b{t) b{t)^ b{t) 



1 



8nG 



-2 



b(t) b{t)^ b(ty^ 



+ A 



+ A 



(7) 
(8) 
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Pit) 



SttG 



a{t) b{t) a{t) bit) 



(9) 



This set of three equations can be solved for a(t), e{t) and P{t) once the initial conditions 
are fixed and the equation of state of the fluid is provided, i.e. P{t) = |e(t) for a radiation 
dominated fluid or P{t) = for a matter dominated fluid. 

In order to get analytical insight on the solutions of this set of equations, let us start 
by considering the limit Tq ^ ^ b(t) in which the comoving radius of the hyper cylinder 
expands/contracts very fast (indeed, much faster than light). In this limit the term r^^/6(t)^ 
proportional to the curvature of the hypercylinder can be neglected and, hence, we can set 
a(t) = bit) and the set of three equations reduces to the usual two equations of the standard 
FLRW flat cosmology. 



e(t) 



SttG 



aiif_ A 
^ 3 



Pit) 



8nG 



a(t) a(t)^ / 
ait) a(t)2 ^ 



In this framework the energy density of the universe evolves roughly as 



(10) 
(11) 



eit) 



,matter 



(12) 



a^it) a%t) ' 

(this expression is only qualitatively true, though, as it does not take into account that at 
high temperatures massive species behave relativistically and contribute to the radiation en- 
ergy density, while at lower temperatures they can decay into lighter modes) and, therefore, 
the Hubble expansion rate Hit) = ait) /ait) evolves as 



H'it) 



SttG / e^ rad , ^0,matter \ ^ 



(13) 



3 \a^it) a%t) J 3' 
Hence, the radiation density controls the expansion rate of the universe at very early times 

and the matter density takes over later on. If the cosmological constant is not zero, A ^ 
0, it eventually becomes dominant at some point. In the concordance ACDM model the 
cosmological constant becomes dominant only at present time. 

Let us now come back to the full set of equations ((7[ [sj [o]) that describe the dynam- 
ics of the hypercylindrical cosmological framework. For the sake of simplicity we rewrite 
the equations in terms of r = Hot, where Hq is the current best observational estimate 
(within the framework of the standard FLRW cosmology) of the Hubble expansion rate at 
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present. As we are interested in exploring the possibility that the novel term proportional 
to the curvature of the hypercylinder r^"^ /h{t)'^ could mimic the contribution of a non-zero 
cosmological constant, from now on we fix A = in these equations: 



Hi 

8nG 



Pi 



1 da{"i 



1 dh{T) 



a{T) dr I \ &(r) dr 



' 1 dh{r 
h{r) dr 



Hi 

SttG 



1 d%{i 



+ 



1 db{T) 



{Ho'/ro 



6(r) rfr2 U(r) dr 



P(t) 



Hi 



1 d^a{T) 1 d'^bir) 



+ 



1 (ia(r) 



1 dbirY 



(14) 



(15) 



(16) 



a(r) dr'^ 6(r) dr^ \'^(t) dr J \fe(r) dr 
For ro the curvature term is negligible at very early times and, therefore, the hy- 

percylindrical cosmological framework is indistinguishable from the standard FLRW flat 



cosmology (10), (11) during the early stages of the history of the universe. Only at present 



time the curvature term becomes non-negligible and the hypercylindrical cosmological model 
diverges from the standard FLRW fiat cosmology. In particular, the two scale factors do 
not remain equal anymore a{t) ^ b{t) and, therefore, neither the Hubble expansion rates 
Ha{t) = d{t)/a{t) and Hi,{t) = b(t)/b(t) remain equal. 

The late evolution of the scale factors within this novel cosmological framework is plotted 
in Fig. 2 (a(t) - red hne; b{t) - blue line) as a function of cosmological time for Hq^ /tq = 0.53 
and initial conditions that correspond to those of the standard FLRW cosmology at the 
instant of decoupling. The solution is compared to the evolution of the scale factor in the 
concordance ACDM model (black line) and in the matter dominated fiat FLRW cosmology 
(grey line) for the same given initial conditions. These initial conditions were normalized such 
that within the framework of the matter dominated fiat FLRW cosmology the scale factor 
at time t = Hq^ would be equal to 1. For the sake of simplicity, in all the three theoretical 
frameworks the contribution of radiation was neglected since the instant of decoupling (that 
is, we set P{t) = 0). It can be seen that the evolution of the scale factor along the axial 
direction in the hypercylindrical model fits very well with the phenomenological ACDM 
model, without any need to postulate any kind of dark energy. The value Hq^/tq = 0.53 
used here is only an arbitrary choice and definitely not the output of any optimization fit. 
Such an optimization will be carried separately. Let us remark that b{t) is the scale factor 
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upon a 2D sphere that contains the plane of our Galaxy, while a{t) describes the dynamics 
along the orthogonal spatial direction. In Fig. 3 is plotted the evolution of the energy 
density in this hypercylindrical model in units of ecr = 

The so-called flatness problem of standard FLRW cosmology gets immediately solved 
within the hypercylindrical cosmological framework that we have just described, for its 
spatial manifold is geometrically flat for any value of the cosmological energy density eo or 
the radius vq [6]. Hence, the model can naturally fit the location of the acoustic peaks in the 
spectrum of primodial cosmic microwave temperature anisotropics over the sky. Moreover, 
the hypercylindrical universe can also accommodate recent claims of a preferred direction in 
the sky map of primordial temperatures anisotropics in the cosmic microwave background 

In Fig. 4 we plot the evolution of the same hypercylindrical model of Fig. 2 and Fig. 
3 well beyond the present time. The scale factor b{t) collapses in finite time and drives 
the exponential expansion of the scale factor along the axial direction a{t). This behavior 
suggests that the collapse of large additional extra dimensions in the very early stages of 
the universe might have driven a kind of inflation along the remaining present three spatial 
dimensions. This issue will be explored separately. 

The hypercylindrical solution can also be smoothly matched to, for example, a closed 
FLRW universe (hypersphere) with comoving radius tq. We can then speculate about about 
the universe as a tree (see Fig. 5) of hierarchical structures. 
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FIG. 1: Comoving spatial sheet associated to metric ([T]). Each horizontal circle represents indeed 
a 3D sphere. 
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FIG. 2: Evolution of the scale factor within three different cosmological models for identical initial 
conditions at the instant of decoupling: a) matter dominated flat FLRW universe, ilm = 1 (gray 
line); b) ACDM with flm = 0.3 and JIa = 0.7 (black line); c) hypercylindrical model with Hq^ /tq = 
0.53, red line - evolution of the scale factor along the axial direction, blue line - evolution of the 
scale factor along the orthogonal 2D manifold. The purple line signals the present age of the 
universe. The hypercylindrical model fits well with the phenomenological ACDM model without 
any need to postulate any kind of dark energy. 
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FIG. 3: Evolution of the energy density in the hypercyhndrical model. The purple line signals the 
present age of the universe. 
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FIG. 4: Evolution of the scale factor well beyond the present time within three different cosmo- 
logical models with identical initial conditions at the instant of decoupling: a) matter dominated 
flat FLRW universe, = I (gray hue); b) ACDM with = 0.3 and Qa = 0.7 (black hue); c) 
hyper cylindrical model with H^^ /r^ = 0.53, red line - evolution of the scale factor along the axial 
direction, blue line - evolution of the scale factor along the orthogonal 2D manifold. The purple 
line signals the present age of the universe. 
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FIG. 5: Hypothetical structure of the universe. Our patch of the universe might be one of the 
hypercyhnders that asymptotically match into a larger similar structure. Smaller hypercylinders 
could also evolve from our visible patch. 



